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The moving puncture method is analyzed for a single, non-spinning black hole. It is shown that 
the puncture region is not resolved by current numerical codes. As a result, the geometry near the 
puncture appears to evolve to an infinitely long cylinder of finite areal radius. The puncture itself 
actually remains at spacelike infinity throughout the evolution. In the limit of infinite resolution the 
data never become stationary. However, at any reasonable finite resolution the grid points closest to 
the puncture are rapidly drawn into the black hole interior by the F-driver shift condition. The data 
can then evolve to a stationary state. These results suggest that the moving puncture technique 
should be viewed as a type of "natural excision" . 



With the moving punctures technique, vacuum black 
holes are modeled on a numerical grid with M 3 topology. 
The interior of each black hole contains an asymptotic re- 
gion that is represented as a point or "puncture" where 
the gravitational field diverges. To be precise, the ini- 
tial data are chosen to be conformally flat, g ao — ip^gab, 
where g a b is the flat metric in Cartesian coordinates. The 
conformal factor ip behaves like 1/r at each puncture 
where r is the coordinate distance from the puncture 
point. The numerical grid is constructed in such a way 
that the punctures do not coincide with any grid point. 
The initial data are evolved with the BSSN formulation 
of the evolution equations [TJ [2] along with "1+log" slic- 
ing, 
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Here, a is the lapse function, j3 a is the shift vector, and 
K is the trace of the extrinsic curvature. The confor- 
mal connection functions are defined by T a = ~dbg ab 
with det(g ao ) — 1. Variants of these gauge conditions 
are obtained by dropping the advection term (3 a d a a in 
Eq. (JT|) or replacing one or more of the time derivatives 
in Eqs. Q by d/dt - /? a <9 Q .[3j [4] 

The moving puncture method was discovered through 
careful numerical experimentation [SJ [3] and has been 
applied successfully by a number of groups. (See, for 
example, Refs. [H [3 E El |10l CH [12l [T3].) Recent 
analytical and numerical work on single, nonspinning 
(Schwarzschild) black holes was carried out with the pur- 
pose of clarifying how the geometry evolves in the neigh- 
borhood of a puncture. |14[ ITS] These works focus on the 
behavior of the numerical slices in a three-dimensional 
evolution code, and find that the puncture evolves from 
a ip ~ 1/r singularity to a l/\/r singularity as the space- 
like slice evolves to a stationary state. Because the phys- 
ical metric behaves like ip 4 (dr 2 + r 2 d£l) with -0 ~ Vv^ 



near the puncture point r = 0, the geometry repre- 
sents a three-dimensional "cylinder" of infinite proper 
length and finite areal radius. Stationary foliations with 
these properties have been described numerically and 

analytically. HI [J51H7] 

The purpose of this report is to extend these insights. 1 

show that the time required for the puncture to evolve to 
a 1/ %fr singularity increases with increasing resolution. 
In the limit of infinite resolution the puncture would re- 
tain its 1/r character throughout any finite-time evolu- 
tion. The behavior of the puncture that is seen in nu- 
merical studies can be misleading because the puncture 
region is not resolved. 

The geometrical picture of puncture evolution that is 
presented here applies to a single, nonrotating black hole. 
It is not clear whether the same description applies to 
black holes with spin and linear momentum. For punc- 
ture initial data the conformal factor diverges like 1/r at 
each puncture. [TS] This suggests that initially, the evolu- 
tion of the geometry near each puncture point will be in- 
dependent of the spin and momentum. As the conformal 
factor changes, the geometry in the vicinity of a puncture 
could develop a dependence on spin and momentum. 

The discussion will be presented using the notation and 
terminology of the Kruskal-Szekeres diagram, shown in 
Fig. 1 . [TU] The Kruskal-Szekeres diagram represents the 
maximal analytic extension of a Schwarzschild black hole. 

The black hole interior is region II, and we live in the 
asymptotically flat exterior region I. Region III is the 
"other" asymptotically flat region, separated from I by 
a wormhole. The heavy curves are the future and past 
singularities. The axes u and v are Kruskal-Szekeres co- 
ordinates in which the metric reads 
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Each point in the u—v plane is a sphere of areal radius R, 
where R is defined by 



2M 



,R/2M 
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The lines at ~ 45° are the horizons with R = 2M. The 
areal radius is greater than 2M in the exterior regions I 
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FIG. 1: Kruskal-Szekeres diagram of a Schwarzschild black 
hole. 



and III. The areal radius is less than 2M in regions II 
and IV. 

The curved arrows in Fig. 1 show the directions of a 
Killing vector field, which is taken as future pointing in 
our region I. Note that any slice can be extended to a 
stationary foliation by Lie transport along the Killing 
vector field. In the Hamiltonian or 3+1 language, any 
spacelike slice can be evolved into a time-independent 
history if the lapse and shift are chosen such that the 
time flow vector field i M = an 11 + /3 M coincides with the 
Killing vector field. (Here, n M is the unit timclike nor- 
mal to the slices and /3 M is the shift vector expressed as 
a spacetime vector field.) Note that we cannot have a 
stationary foliation with a > everywhere if the slices 
extend into both regions I and III. The foliations found 
in Refs. [TUEOH?] are examples of stationary foliations 
that begin in region II and end in region I, extending 
from u ~ — oo, v ~ co to u ~ co, v ~ finite. For any 
such foliation whose slices asymptotically approach an 
R = const curve in region II (a hyperbola in Kruskal- 
Szekeres coordinates), the slice geometries will approach 
that of an infinitely long three-dimensional cylinder. The 
stationary foliation found in Refs. [TBI EZ] satisfies the 
maximality condition K = 0, while the stationary folia- 
tion discussed in Ref. [H] satisfies the stationary "1+log" 
slicing condition 2aK — (3 a d a a. 

The initial data for puncture evolution of a 
Schwarzschild black hole consist of the spatial metric in 
isotropic coordinates, 

ds 2 = (1 + M/(2r)) 4 (dr 2 + r 2 dtt 2 ) , (5) 

and vanishing extrinsic curvature. These are the data for 
the v — slice of the Kruskal-Szekeres diagram. For the 
initial data the relationship between u and the isotropic 
radial coordinate r is given by 




The positive u axis is covered by M/2 < r < co and the 
negative u axis is covered by < r < M/2. Initially the 
puncture point r = coincides with spacelike infinity in 
region III, u — -co. 

The key question is this. Can the initial data v = 
evolve to a slice that begins in region II and ends in re- 
gion I? The answer is no. We can see this most clearly 
by keeping in mind that, in Hamiltonian language, the 
momentum constraint is the generator of spatial diffco- 
morphisms. As a consequence the time derivatives in the 
3+1 evolution equations appear only in the combination 
d± = d/dt — Cp, where Cp is the Lie derivative along the 
shift vector f3 a . Let us assume that the slicing condition 
has this same property: time derivatives appear only in 
the combination d±. When this is the case, the evolu- 
tion equations and slicing condition together determine 
the history of the (intrinsic and extrinsic) geometry of 
space. The lapse function specifies the point wise proper 
time separation between instances of space (or "slices" ) 
at successive times. The choice of shift vector does not 
affect the spatial geometry. The role of the shift vector 
is to determine how the spatial coordinates are carried 
from one slice to the next. In a finite difference calcula- 
tion the shift vector determines how the distribution of 
grid points is carried from one slice to the next. 

Consider the evolution of the initial geometry repre- 
sented by the v = slice using 1+log slicing, Eq. ([I]). 
We can ignore how the grid points are distributed ini- 
tially, and how they are redistributed in time by the 
shift vector. Now observe that the maximally extended 
Schwarzschild black hole is symmetric under reflections 
about the Kruskal-Szekeres v axis. Thus, the initial met- 
ric expressed in Kruskal-Szekeres coordinates (the metric 
of Eq. ([3) with v = 0) is invariant under u — > —u. The 
initial metric expressed in isotropic coordinates Eq. ^ 
is invariant under r — > M 2 /(4r). This reflection symme- 
try implies that the initial spatial geometry in region III 
is identical to the initial spatial geometry in region I. As 
long as the initial choice of lapse function is also reflection 
symmetric, then the geometry will remain symmetric as 
it evolves in time. If the part of the slice that extends to 
spatial infinity in region I does not evolve into region II, 
then the part of the slice that extends to spatial infinity 
in region III also does not evolve into region II. 

In fact, neither end of the initial data slice v — can 
evolve entirely into the black hole interior (region II) in 
finite time. To see why this is the case, consider a se- 
quence of spacelike slices whose members are separated 
by finite proper time, and whose initial slice is v = 0. 
Starting from a fixed point in region III on the initial 
slice, we can define a timelike trajectory of finite proper 
length that is orthogonal to each of the slices. Now take 
the limit as the fixed initial point moves to spacelikc infin- 
ity. A timelike trajectory cannot extend from spacelikc 
infinity to the black hole interior in finite proper time. 
This is clear from the Penrose diagram [19] , which shows 
that a curve connecting spacelike infinity to the black 
hole interior must be null or somewhere spacelike. It fol- 
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lows that the timelike trajectory remains in region III. 
The spacelike slices must also extend into region III. 

The argument above shows that the puncture point, 
which begins at spacelike infinity in region III, must re- 
main in region III throughout its evolution. All of the 
slices evolved from v = begin in region III and end in 
region I. With a positive lapse these slices cannot reach 
a stationary state. 

These results are confirmed by numerical tests. The 
tests described below are carried out with a BSSN code 
that assumes spherical symmetry. The code is based 
on a generalization of the BSSN formulation that does 
not assume det(<? a fc) = l.[20] The evolved data include 
the components g rr and gee of the conformal metric, the 
"chi version" of the conformal factor \ = l/V^i the trace 
of the extrinsic curvature K, the component A rr of the 
trace-free part of the extrinsic curvature, and the con- 
formal connection function r r . The lapse function a is 
determined by the 1+log slicing condition Eq. (JlJ. Each 
of these variables, along with the shift vector component 
(3 r , is a function of time t and radial coordinate r. The 
initial metric is given in Eq. ([5|. The finite difference 
grid is cell centered with grid points r,j = (j — 1 /2) Ar for 
j = 1,2,.... The puncture is located at r = 0. The code 
also tracks the Kruskal-Szekeres coordinates u and v of 
each grid point; this allows the slices to be displayed in 
a Kruskal-Szekeres or Penrose diagram. Further details 
of the 1-D code are described in Ref. [2"T] . 

Figure 2 shows the slice at time t — 3Af as a curve 
just below the future singularity on the Kruskal-Szekeres 
diagram. The slice begins in region III and ends in region 
I. This slice was obtained by evolving the initial data v = 
with initial lapse function a = 1. For this simulation 
the shift vector was set to zero for all time. The reflection 
symmetry about the w-axis is evident. The heavy dots in 
Fig. 2 show the locations of grid points obtained from a 
second simulation, identical to the first except for the use 
of the T-driver shift condition. For clarity of presentation 
only every other grid point is displayed, beginning with 
grid point j — 2. The resolution used for both of these 
runs was Ar = M/50; this is comparable to the highest 
resolutions used in current 3D codes. 

With the T-driver shift, the left-hand side of the 
Kruskal-Szekeres diagram is almost completely devoid 
of grid points by t = 3Af. The asymmetry in the distri- 
bution of grid points occurs for two reasons. First, the 
grid points in the initial data are distributed uniformly 
in r, not in u. With an outer boundary at, say, 100A/, 
only 0.5% of the grid points start in the left half of the 
Kruskal-Szekeres diagram. Second, with the T-driver 
condition the shift vector quickly develops (3 r ~ r behav- 
ior near the puncture. Since the conformal factor diverges 
like 1/r, the magnitude of the shift vector behaves like 
y/{3 r ip 4 g rr l3 r ~ 1/r. Due to its large magnitude at small 
r, the shift vector quickly drives the grid points in region 
III toward the black hole interior. 

The spherically symmetric code differs from current 
three-dimensional black hole codes in one important re- 
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FIG. 2: Slice at t = 3M for vanishing shift (solid curve) and 
F-driver shift (dots). 



spect: with the spherically symmetric code the puncture 
point r = is a boundary of the computational domain. 
In most contexts, the choice of boundary conditions will 
affect the evolution of the fields in the bulk. For puncture 
evolution, however, the results are found to be insensi- 
tive to the choice of boundary conditions at r = 0. This 
is most likely due to that fact that, with the T-driver 
shift condition, the grid points near the puncture quickly 
acquire superluminal speeds as they are drawn into the 
black hole interior. The grid points continue to evolve 
along spacelike trajectories that approximate the orbits 
of the Killing vector field. It follows that at the bound- 
ary r — all of the physical characteristics are directed 
from interior to exterior (toward the left in the Kruskal- 
Szekeres diagram). As a result one would not expect 
conditions imposed at r — to influence the evolution of 
the geometry. 

With the spherically symmetric code we can easily in- 
crease the resolution well beyond the value used in the 
simulations of Fig. 2. This is an effective way to coun- 
teract the sparsity of grid points in region III in the ini- 
tial data. For example, with Ar = M/50 and the outer 
boundary at 100M there are only 25 grid points with 
initial values u < 0. We can increase this by a factor 
of ten or more and still maintain reasonably short run 
times. However, a straightforward increase in resolution 
is not very effective in counteracting the movement of 
grid points caused by the T-driver shift condition. At 
resolutions as high as Ar = A//6400 the shift vector still 
drives all of the grid points into the black hole interior 
within a time of a few M. See Fig. 3. After all of the grid 
points have left region III, the lapse and shift can settle 
into a stationary state where the time flow vector field 
f = an^ + /3 M coincides with the Killing vector field. 

For higher resolutions the time required for all of the 
grid points from region III to be driven into the black hole 
interior increases, as seen in Fig. 3. In the limit of infi- 
nite resolution, the grid points at any finite time would 
stretch from region III to region I, crossing the black hole 
interior like the curve in Fig. 2. At infinite resolution the 



4 




1 1.5 2 2.5 3 3.5 4 

log(M/Ar) 



FIG. 3: The vertical scale is the time required (in units of 
M) for all of the grid points that begin in region III to evolve 
into the black hole interior, region II, with the T-driver shift 
condition. The horizontal scale is the common logarithm of 
M/Ar. The resolutions used for this graph differ by powers 
of 2, ranging from Ar = M/12.5 to Ar = Af/6400. 
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FIG. 4: Logarithm of the lapse a as a function of areal radius 
R (in units of M) at time t — 30M for vanishing shift (solid 
curve) and T-driver shift (dots). 

slices never reach a stationary state. However, the time 
T required for the region III grid points to evolve into 
region II increases very slowly with resolution, roughly 
like T/M ~ 0.58 log(M/Ar) + 1.5. We can compare this 
to the time scale for the lapse and shift to settle to a sta- 
tionary state — at resolutions typical of current codes this 
is about 50M. According to the heuristic expression for 



T/M, an extraordinarily high resolution of Ar ~ M/10 83 
would be required for grid points to remain in region III 
beyond t ~ 50M. 

One of the important properties of the 1+log slicing 
condition is that the lapse function is driven to zero in 
the black hole interior, so the slices naturally avoid the 
physical singularity. Figure 4 shows the common loga- 
rithm of the lapse at time t = 30M plotted as a function 
of areal radius for both u > (right half of the graph) 
and u < (left half of the graph). The continuous curve 
was obtained from a simulation with vanishing shift vec- 
tor and resolution Ar = M/800. The dots were obtained 
from a simulation with the T-driver shift condition and 
resolution Ar = M/50. (For clarity, only every other 
grid point is displayed.) Again we see that the evolu- 
tion of the complete slice is symmetric under reflections 
u — > —u, but for the simulation with T-driver shift the 
distribution of grid points is highly asymmetric. 

With the moving puncture method the region inside 
the black hole horizon is not completely resolved. As a 
matter of principle, the puncture point remains at space- 
like infinity in region III throughout its evolution. How- 
ever, with the T-driver condition the shift vector quickly 
drives all of the region III grid points into the black hole 
interior. The resulting slices are incomplete; the grid 
points terminate in region II due to lack of resolution. 
The lapse function and shift vector adjust to bring the 
data to a stationary state. The stationary state is a por- 
tion of a slice that extends from u ~ — oo, v ~ oo in 
region II to u ~ oo, v ~ finite in region I. 

The lack of resolution that characterizes the puncture 
method is not likely to cause problems for finite difference 
codes at any reasonable resolution. In fact, the results 
here suggest that puncture evolution should be viewed as 
a type of "natural excision" . With black hole excision one 
intentionally removes from the computational domain a 
set of grid points inside the horizon (see, for example, 
Ref. [22] )• In the case of puncture evolution the removal 
of grid points occurs naturally. The advantage over the 
traditional excision method is that puncture evolution 
does not require the introduction of an interior boundary 
where boundary conditions and special finite difference 
stencils must be applied. 
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